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INTRODUCTION 


The  reduced  cohesion  of  grain  boundaries  (GB)  is  often  the  controlling  factor  limiting 
the  ductility  of  high  strength  metallic  alloys.'  This  is  particularly  .so  when  there  is  an  environ¬ 
mental  interaction,  as  in  the  case  of  intergranular  hydrogen  stress  corrosion  cracking. “  Typi¬ 
cally,  hydrogen  drastically  degrades  mechanical  properties  of  metal  alloys,  making  them 
unreliable  and  poses  a  significant  technological  problem. 

Intergranular  embrittlement  in  metals  is  usually  associated  with  prior  segiCgation  of  impurities 
toward  the  GBs.^'^  Impurities  pre.sent  in  bulk  concentrations  of  lO'^-lO”*  atomic  percent  can  result 
in  a  dramatic  decrease  in  plasticity.  The  hydrogen  concentration  is  also  extremely  low;  its  solid 
.solubility  in  bcc  Fe  below  room  temperature  is  only  1  ppm  (lO"^  atomic  percent).^ 

Previous  studies  attempted  to  explain  GB  decohesion  on  the  electron-atom  level  due  to 
hydrogen  by  using  small  atomic  clusters  as  a  simple  model  of  GB.^  Decohesive  effects  of 
hydrogen  were  akso  succe.ssfully  investigated  using  semiempcrical  calculations  and  modelling'^ 
based  upon  the  Embedded  Atom  Method  (EAM). It  was  shown  that  a  single  hydrogen 
atom  in  a  unit  cell  at  a  GB  of  fee  Ni  can  weaken  the  metallic  bonds  across  the  GB  lowering 
the  fracture  stress  by  15%. 

Progre.ss  in  developing  efficient  first-principles  methods  in  the  recent  decade  has  allowed 
more  detailed  treatments  of  a  GB  with  an  impurity.  Calculations  making  use  of  both  two- 
dimensional"’*^  and  supercell’’^  models  have  provided  a  deeper  insight  into  the  mechanisms 
of  impurity  decohesion. 

The  effect  of  hydrogen  on  the  electronic  structure  of  transition  metals  has  also  been  a 
focus  of  numerous  research  papers*"^  mainly  directed  at  an  investigation  of  stable  hydrides. 
However,  since  Fe  does  not  form  a  compound  with  hydrogen,  neither  hydrogen  in  the  bulk 
Fe  nor  in  Fe  GB  has  been  studied  in  detail  thus  far. 


The  present  paper  reports  preliminary  results  of  the  self-consistent,  spin-polarized  calcula¬ 
tions  of  the  electronic  structure  of  H  in  an  atomic  environment  typical  of  a  GB  in  Fe.  A 
more  detailed  analysis  will  be  reported  at  a  later  date. 

Atomistic  relaxation  studies***  have  shown  that  an  interstitial  impurity  atom,  such  a  boron, 
is  likely  to  occupy  the  center  of  a  capped  trigonal  prism  formed  by  Fe  atoms  in  the  core  of 
the  GB.  Hydrogen,  due  to  its  small  atomic  size,  may  po.ssibly  occupy  tetrahedral  voids  in  the 
GB  plane  rather  than  centers  of  trigonal  prisms.  In  the  present  study,  however,  H  has  been 
treated  in  the  same  way  as  B.  C.  P,  and  S  atoms  in  the  Fe  GB.*^ 


A  26-atom  supercell  was  used  as  a  model  of  the  GB  (.see  Figure  1)  with  the  correspond- 
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Figure  1 .  26-atom  supercell  model  of  GB.  a)  The 
hexagonal  supercell;  only  parts  of  the  unit  cell 
adjacent  to  the  GB  are  shown,  b)  The  trigonal 
prism  environment  on  the  GB  planes,  c)  View 
along  [1 12]  direction  of  the  periodic  (111)  plane 
array;  small  and  large  circles  are  Fe  atoms  in  alter 
nating  (110)  planes.  The  capped  trigonal  prisms 
on  the  GB  are  shown.  The  broken  line  rectangles 
centered  at  Fe-7  atoms  show  the  projections  of 
bcc  unit  cells. 


The  two  stacking  faults,  ...BAB...  and  ...CAC...  imitate  two  (111)  23  tilt  GB.  One  of  the 
GBs  is  occupied  by  an  H  atom,  while  the  other  is  "clean;"  an  empty  sphere  (ES)  is  placed 

inside  the  trigonal  prism  of  Fe  atoms  (see  filled  square  in  Figure  1).  There  arc  11  planes  of 

Fe  atoms  between  the  two  GB  planes  (A).  It  is  important  that  information  on  both  the 

clean  GB  and  the  GB  with  H  can  be  obtained  within  the  same  calculation.  In  order  to  find 

the  energy  difference  resulting  from  placing  an  H  atom  on  the  GB.  an  independent  reference 
calculation  on  a  superccll  with  two  clean  GBs  was  also  performed. 

To  make  the  GB  model  more  realistic,  both  the  clean  GB  and  that  with  H  were  "ri'lnv^ri" 
by  minimi^ii.g  energy  ar,  a  fnn^tinii  .T  inierplanar  spacings.  I  he  Finnir  i,mciaii  ^  version  ot 
EAM  was  used;  the  energy  of  an  H  atom  was  found  as  a  function  of  local  electron  charge 


density  at  the  H-site  due  to  surrounding  Fe  atoms  from  the  so-called  Effective  Medium  The¬ 
ory  (EMT).‘^ 

Spin-polarized  scalar-relativistic  calculations  were  done  using  the  Linear  Muffin-Tin  Orbitals 
(LMTO)  method. The  exchange-correlation  functional  of  von  Barth  and  Hcdin"®  and  ihe 
frozen  core  approximation  were  used. 

RESULTS  AND  DISCUSSION 

Figure  2  shows  the  relaxed  interplanar  distances  for  both  the  clean  GB  and  the  GB  with 
H.  An  interesting  feature  of  the  clean  GB  is  that  the  distance  between  the  second  and  third 
planes  is  only  half  of  the  (111)  interplanar  distance  in  bulk  bcc  Fe.  The  consequences  of 
this  feature  will  be  discussed  below. 


Figure  2.  "Relaxed"  distances  between  (1 1 1)  GB  planes  num¬ 
ber  i  and  1-1  for  the  clean  GB,  and  GB  with  hydrogen  as  a  func¬ 
tion  of  the  plane  number  I.  The  GB  plane  is  plane  1 ;  atoms  Fe-1 
belong  to  plane  1 ,  Fe-2  to  plane  2,  etc.. 


Figures  and  4  show  the  electronic  densities  of  stales  (DOS)  in  units  of  stalcs/Ry  versus 
energy  in  mRy.  The  zero  of  energy  corresponds  to  the  Fermi  energy  (Ej.)  of  perfect  bcc  Fe. 
The  arrows  show  the  Eyr  for  the  clean  GB  and  GB  with  H. 

In  Figure  3  we  compare  the  DOS  of  a  bulk  bcc  Fe  atom  with  the  site-projected  DOS  of 
atom  Fe-7  on  the  central  plane  between  the  two  GBs.  The  two  DOS  plots  are  very  much 
alike,  which  suggests  that  the  supercell  size  is  sufficient  for  minimal  interaction  of  the  two 
GBs.  It  is  observed,  however,  that  the  d-bands  of  the  Fe-7  atom  are  somewhat  narrowed, 
sinc’f"  fbc  distances  between  Fe-7  and  its  neighbors  are  slightly  increased  due  to  the  "tail"  of 
the  elastic  "relaxation  wave;"  tiic  latter  "dumps  out"  only  after  13  to  15  crystal  planes  from 
the  GB  (sec  Figure  2). 
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Figure  3.  DOS  for  the  bulk  bcc  Fe  atom  (upper  plot)  and  Fe-7  (center) 
atom  in  the  supercell  (lower  plot). 


Figures  4a  through  4c  compare  the  site-projected  DOS  for  atoms  Fe-l,  Fe-2,  and  Fc-3 
for  the  clean  GB  (upper  plots)  and  the  GB  with  H  (lower  plot).  The  bottom  ("negative") 
parts  of  the  lower  plots  show  the  site-projected  DOS  within  the  Wigncr-Seitz  (WS)  sphere  of 
the  H  atom.  These  DOS  are  identical  for  the  three  plots. 

The  DOS  of  the  Fc-1  atom  on  the  clean  GB  is  strongly  distorted  compared  to  that  for  a 
bulk  atom.  This  can  be  easily  understood  since  relaxation  of  the  clean  GB  brings  the  (111) 
planes  of  the  Fc-2  and  Fe-.3  atoms  very  close  to  each  other  as  we  mentioned  above  (the  intcr- 
pianar  distance  is  only  0.407A.  compared  to  0.825A  in  the  bulk).  It  is  only  the  interaction  of 
the  Fc-3  plane  with  the  outer  planes  Fc-4,  Fc-5,  and  Fc-6  which  prevents  it  from  collapsing 
into  an  a>-phasc  configuration.  In  fact,  the  DOS  for  the  Fe-l  atom  on  the  clean  GB.  as 
shown  in  Figure  4a,  is  very  much  like  the  DOS  for  the  cy-phase.^'  Off  the  GB  plane,  in  Fe-2 
and  Fe-3,  the  DOS  acquire  more  bcc-like  features. 

The  presence  of  H  on  the  GB  pushes  the  Fe-2  and  Fe-3  planes  apart  making  the  Fe-l 
DOS  more  bcc-like.  In  Fc-2  and  Fe-3  the  bcc  shape  is  further  restored.  The  d-band  width 
behaves  as  one  could  have  expected;  i.c.,  following  the  lattice  distortion  pattern,  being  in  gen¬ 
eral,  narrower  on  the  GB  with  H.  A  number  of  peaks  around  Ep  on  both  the  GBs  are  char¬ 
acteristic  of  localized  GB  states. 
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Figure  4  a-c.  The  site-projected  DOS  for  atoms  Fe-1 ,  Fe-2, 
and  Fe-3  for  the  clean  GB  (upper  plots)  and  the  GB  with  H 
(lower  plots).  The  bottom  ("negative")  parts  of  the  lower 
plots  show  the  site-projected  DOS  within  the  WS  sphere  of 
the  H  atom  (same  for  the  three  plots). 
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When  H  is  placed  on  the  iron  GB,  Fe-H  chemical  bonds  are  formed,  as  in  the  case  of 
other  impurities.'^  These  are  the  consequence  of  electron  hybridization.  Since  the 
ground  state  of  the  electron  in  an  isolated  H  atom  (frj  in  Figures  4a  through  4c)  corre¬ 
sponds  to  an  energy  range  where  the  Fe  pd-components  are  extremely  weaK,  only  a  weak  s- 
pd  hybridization  can  be  expected.  The  latter  results  in  forming  a  bonding  band,  with  states 
of  mostly  s-type  within  the  WS  sphere  of  H,  and  of  mostly  pd-type  in  the  WS  spheres  of  the 
Fe  atoms.*  A  simple  estimate  shows  that  an  antibonding  component  of  the  hybrid  should 
appear  in  the  energy  range  about  (-160  to  -180)  mRy.  In  fact,  one  can  see  a  series  of  tiny 
peaks  superimposed  on  the  "bcc-like"  dip  in  the  DOS  plots  in  this  energy  range.  Unlike  a 
more  typical  case  when  impurity  antibonding  states  are  located  above  Ep,  these  antibonding 
states  are  filled,  thus  resulting  in  a  decrease  of  stability. 

The  strongest  bonding  d-component  is  localized  on  the  Fe-1  atoms  even  though  the  Fe-3 
atom  is  closer  to  the  H  site;  the  numbers  of  d-elcctrons  in  these  states  are  0.069  and  0.045, 
respectively.  This  indicates  that  the  s-d  hybrid’s  a-component  spreading  out  in  the  GB  plane 
is  stronger  that  a  <5-component,  penetrating  across  the  GB. 

A  question  of  importance  is  what  happens  with  the  H  electron?  It  has  been  generally 
believed^’'’'’^^  that  in  group  7  and  8  transition  metals  H  donates  its  electron  to  a  host’s  valence 
d-band.  The  results  of  Reference  13  for  the  Fe  GB  with  B,  C,  P,  and  S  showed  there  was  a 
strong  charge  redistribution  among  the  atomic  spheres,  but  the  present  results  show  that  this  is 
not  the  case  for  H.  Table  1  presents  the  *otal  numbers  of  electrons,  n,  within  atomic  WS 
spheres  (.Sq  are  the  sphere  radii  in  a.u.),  their  s-,  p-,  and  d-components,  and  the  magnetic 
moments,  m  (in  //B/atom).  Though  the  concept  of  charge  transfer  is  quite  ambiguous,  neverthe¬ 
less,  some  of  the  features  of  the  charge  redistribution  pattern  shall  be  discussed. 

Since  the  WS  radii  of  the  Fe  atoms  of  "relaxed"  GB  with  and  without  H  are  different,  it 
is  difficult  to  compare  the  contents  of  s-  and  p-  electrons  in  the  corresponding  WS  spheres. 
Comparison  of  the  d-electron  numbers  is  more  meaningful,  since  the  d-slatcs  are  more  local¬ 
ized,  and  the  d-electron  content  should  not  be  sensitive  to  the  WS  radii.  It  can  be  seen 
from  Table  1  that  upon  placing  H  on  the  GB  the  number  of  electrons  localized  on  Fe-1,  Fe- 
2,  and  Fe-3  atoms  increase,  respectively,  by  0.186,  0.027,  and  0.010  electrons/atom.  The  calcu- 
laUons  show  that  the  numbers  of  d-electrons  localized  in  the  impurity  bonding  bands  are 
0.069,  0.030,  and  0.045,  respectively;  then,  the  numbers  of  d-electrons  in  the  "valence"  part  of 
the  spectrum  changes  by  +0.117,  -0.003,  and  -0.035.  Thus,  the  total  change  in  the  d-electron 
content  per  H  atom  (based  upon  lFe-1,  2Fe-2,  and  2Fe-3  atoms  per  each  H  atom  on  the 
GB,  as  shown  in  Figure  1)  is  0.260,  of  which  0.217  d-electrons  arc  localized  in  the  Fe  bond¬ 
ing  bands,  while  the  valence  d-band  population  increases  by  only  0.041  electron.  At  the  same 
time,  there  is  a  "deficiency"  of  0.238  electrons  in  the  WS  sphere  of  the  H  atom. 


*A.s  fix’  the  ^components  of  the  t'e  impunty  bands,  they  are  merely  "tails”  of  the  s-cicetron  dlsiributkin  in  the  llydmgcn  WS  sphere. 
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Table  1. 


Atom 

So 

ns 

np 

na 

n 

m 

■ 

0.753 

0.045 

0  004 

0  000 

0.049 

0  002 

Fe-1 

2  834 

0  727 

0935 

6.689 

8.350 

2  595 

Fe-2 

2.519 

0.593 

0  638 

6  490 

7722 

2.043 

Fe-3 

2.654 

0625 

0.732 

6572 

7  929 

2  339 

Fe-4 

2  790 

0721 

0919 

6.694 

8.335 

2.387 

Fe-7 

2.661 

0.615 

0  730 

6  505 

7  849 

2  340 

bcc,  bulk 

2.661 

0.641 

0.784 

6.575 

8.000 

2  298 

Fe-4 

2.778 

0.710 

0.888 

6  664 

8  262 

2.395 

Fe-3 

2.617 

0.614 

0.710 

6  582 

7  907 

2.097 

Fe-2 

2  588 

0.620 

0.649 

6518 

7.788 

2232 

Fe-1 

2  961 

0.761 

0.992 

6  875 

8  628 

2.501 

H 

1.326 

0.697 

0.059 

0  006 

0762 

-0017 

In  spite  of  the  ambiguity  of  charge  transfer,  the  atomic  sphere  approximation  (ASA)  used 
in  the  LMTO  method  enables  one  to  draw  meaningful  conclusions  regarding  electrostatic  inter¬ 
actions  at  the  GB.  Our  calculations  show  (see  Table  1)  that  the  change  in  net  snhere 
charges  upon  placing  the  H  atom  on  the  GB  changes  the  electrostatic  interaction  between  the 
spheres.  On  the  clean  GB,  the  electronic  charge  localized  within  the  ES  has  an  electrostatic 
repulsion  with  the  excess  charge  on  Fe-1  and  attraction  with  the  charge  deficient  Fe-2  and  Fe-3. 
Thus,  the  Fe-3  atoms  across  the  GB  are  pulled  towards  each  other  by  their  attraction  to  the 
ES.  On  the  other  hand,  on  the  GB  with  H,  due  to  the  electronic  deficiency  in  the  H 
sphere,  there  is  an  attraction  between  H  and  Fe-1,  and  repulsion  between  H  and  Fe-2  and  Fc-3. 
Thus,  the  Fe-3  atoms  across  the  GB  undergo  repulsion  from  each  other. 

Finally,  a  few  words  about  the  effect  of  H  on  the  ferromagnetic  moment,  m  (in  ,«jj/atom).  of 
the  nearby  Fe  atoms.  First  oi  all,  the  enhancement  of  m  (2.595)  of  the  Fe-1  atom  on  the 
clean  GB  plane  with  respect  to  that  of  the  center  Fe-7  atom  (2.340),  or  the  bcc  bulk  atom 
(2.298),  is  similar  to  that  expected  around  a  two-dimensional  defect.  Enhancement  of  m  on 
Fe  free  surfaces  has  been  consistently  predicted. Similar  to  the  case  of  the  free  surface, 
the  origin  of  the  moment  enhancement  in  GB  is  the  emergence  of  the  localized  GB  states  in 
the  d-component  of  the  DOS  at  Ep.  This  enhancement  is  slightly  weakened  (2.501)  upon  add¬ 
ing  an  H  atom  to  the  GB,  possibly  due  to  the  narrowing  of  the  d-band  and  thus  decreasing 
the  d-DOS  at  Ep.  A  more  considerable  suppression  of  m  is  found  for  the  Fe-3  atom  (2.097 
versus  2.339  on  the  clean  GB).  The  moment  on  the  Fe-2  atoms  is  only  slightly  enhanced. 

As  for  the  H  atom,  it  appears  to  have  a  very  low  m  (0.0173),  due  to  spin  pt)iarization 
of  s-electrons,  antiferromagnetically  coupled  to  the  Fe  magnetic  moments. 
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CONCLUSIONS 


It  has  been  found  that  in  a  typical  Fe  GB  environment,  hydrogen  lias  a  very  weak  bond¬ 
ing  with  the  iron  atoms  both  in  th  i  GB  plane  and  across  the  GB;  the  latter  being  relatively 
weaker  than  the  former.  A  weak  electrostatic  interaction  and  a  suppression  of  ferromagnetic 
moments  of  the  iron  atoms  across  the  GB  also  contribute  to  decohesion.  The  analysis  of  the 
d-charge  redistribution  showed  that,  contrai^  to  a  general  belief,  hydrogen  on  the  GB  gives  a 
negligible  contribution  of  electron  charge  to  the  i^on  valence  d-bands.  These  results  should 
be  important  for  understanding  the  deenhesion  mechanisms  which  can  produce  hydrogen 
embrittlement  in  iron  and  iron-base  alloy.; 
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